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Abstract 

In this paper, we present the results of Monte Carlo simulations for two popular techniques of long- 
range correlations detection - classical and modified rescaled range analyses. A focus is put on an 
effect of different distributional properties on an ability of the methods to efficiently distinguish 
between short and long-term memory. To do so, we analyze the behavior of the estimators for 
independent, short-range dependent, and long-range dependent processes with innovations from 
8 different distributions. We find that apart from a combination of very high levels of kurtosis 
and skewness, both estimators are quite robust to distributional properties. Importantly, we show 
that R/S is biased upwards (yet not strongly) for short-range dependent processes, while M-R/S 
is strongly biased downwards for long-range dependent processes regardless of the distribution of 
innovations. 

Keywords: rescaled range analysis, modified rescaled range analysis, Hurst exponent, long-term 
memory, short-term memory 
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1. Introduction 

Examination of long-range dependence has become more frequent in various fields in recent 
years. The portfolio of disciplines has grown wide and varies across finance [12, 3, 19], cardiology 
[39, 34], vascular science [21], DNA sequences [18], pollution examination [38], river flow fluctu- 
ations [30], geochemistry [41], earthquakes waiting times [9], climatology [33] and geomagnetism 
[16]. Hurst exponent H, < H < 1, is a measure of long-range dependence for a stationary pro- 
cess with an autocorrelation function p{k) decaying as p{k) oc k 2H ~ 2 with a lag k —> oo. H = 0.5 
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indicates two possible processes - either an uncorrclated or a short-range dependent process. Au- 
tocovariances of the uncorrelated process are insignificant for all non-zero lags whereas for the 
short-range dependent process (e.g. ARIMA processes, for details, see [8]), autocovariances are 
significant at low lags and insignificant at high lags and decay exponentially. If H > 0.5, the 
process has positive correlations at all lags and is said to be long-range dependent with positive 
correlations or persistent. On the other hand, if H < 0.5, the process is long-range dependent 
with negative correlations or anti-persistent. Autocovariances of long-range dependent processes 
are hyperbolically decaying and are either non-summable for the persistent or summable for the 
anti-persistent processes [7, 22, 14, 26, 5]. 

In a majority of applied research papers focusing on long-range dependence in a time series, 
the results are interpreted on a basis of simple comparison of the estimated Hurst exponent H 
with the value of 0.5. If the estimate deviates from 0.5, the studied process is claimed to be 
long-range dependent. However, the estimators of Hurst exponent are usually either biased or 
have high variance for finite samples as shown in several studies (e.g. [37, 11, 15, 23, 6, 20]). The 
estimators can also be biased by a presence of short-range dependence in an underlying process 
[29, 23]. However, distinguishing between the short and the long-range dependence is crucial for 
time series analysis as each type of dependence implies different features of the scries (see [7] for 
discussion). To add to the discussion, we test two methods which are closely related and one of 
them was constructed to be robust against the short-range dependence - the rescalcd range analysis 
of Hurst [17], which is also called the "classical", and the modified rescalcd range analysis of Lo 
[23]. Finite sample properties of the two methods and efficiency of the distinguishing between 
the types of dependence have been partially discussed by Teverosky et al. [35]. However, the 
simulated processes always have normally distributed innovations. Moreover, Teverosky et al. [35] 
tested the ability of the modified rescaled range analysis to adjust the rescaled ranges from the 
classical form of the method rather than its ability to estimate Hurst exponent H . To fill the 
gap, we test the robustness of the rescaled range analyses for independent, short-range dependent 
and long-range dependent processes with innovations from 8 different distributions - standard 
normal, log-normal, Cauchy, log-i, gamma, inverse gamma, Laplace and log-Laplace. By a choice 
of distribution, we discuss on the robustness of both methods to skewness, kurtosis and infinite 
moments as well as to different types of memory. 

The paper is structured as follows. In Section 2, we present and describe both techniques in 
detail. In Section 3, we present important parameters of Monte Carlo simulations methodology 
and describe the chosen distributions. In Section 4, we show results of the Monte Carlo simulations 
for time series lengths from 512 to 16,384 observations for both tested methods and different types 
of distributions of the underlying process. We find four important outcomes - extreme excess 
kurtosis combined with extreme skewness causes the estimates to be downward biased; M-R/S 
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in its simplest way is not able to distinguish between short and long memory efficiently; R/S is 
biased by the short-range dependence in the series but in a lesser way than standardly claimed; 
and apart from very extreme cases of high skewncss and lcptokurtosis, both methods are quite 
robust to different distributional properties. Section 5 concludes. 

2. Classical and modified rescaled range analysis 

The classical rescaled range analysis (R/S) was developed by Edwin Hurst while working as an 
engineer in Egypt [17] and was later introduced to the financial time series by Benoit Mandelbrot 
[26, 27, 25]. For time series X t of length T with t = 0, 1, . . . ,T, Hurst exponent H can be estimated 
from behavior of rescaled ranges of the series. In the procedure, one takes the increments of the 
time series xt = Xt — Xt-i with t = 1, . . . ,T and divides them into N adjacent sub-periods 7 n , 
for n = 1, . . . , N, of length v while Nv = T. For each sub-period, one calculates an average value 
and constructs a new series of accumulated deviations from the arithmetic mean values - a profile 
y t - The procedure follows in a calculation of the range of the profile Rj n = maxj n (y t ) — minj n (y t ) 
and a standard deviation S7„ of the original increment time series for each sub-period. Each range 
is standardized by the corresponding standard deviation and forms the rescaled range so that the 
average rescaled range for the given sub-period length (R/S) v is calculated [32]. Rescaled ranges 
scale as (R/S) v oc v H . 

R/S was shown to be sensitive to hctcroskcdasticity and short-range dependence in the un- 
derlying process [24, 1]. To deal with the problem, Lo [23] proposed the modified rescaled range 
analysis (M-R/S). The procedure differs in the definition of the standard deviation which deals 
with both heteroskedasticity and short-range dependence. The modified standard deviation is 
defined with a use of auto-covariance "fj of the selected sub-interval I n for lag j = 1, . . . , £ as 



Thus, R/S turns into a special case of M-R/S with £ = 0. The most problematic and also the 
crucial issue of the new standard deviation measure is the number of lags £ which are used for its 
estimation. On one hand, if the chosen lag is too low, it omits the lags which may be significant 
and therefore the estimates of the rescaled ranges and Hurst exponent are still biased. On the 
other hand, if the used lag is too high, the estimates of the rescaled range differ significantly from 
the true values [35, 36]. 

One can either estimate the rescaled ranges and Hurst exponents for arbitrary values of £ [40, 2] 
or use an automatic estimator of £ for each sub- interval of the length v. Lo [23] suggests to set the 
optimal lag £* based on the sample first-order autocorrelation coefficient p(l) of the sub- interval 
/„ as (where [J is the nearest lower integer operator) 




(1) 
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Of course, Eq. 2 assumes that the underlying process is some specification including AR(1) 
process. Teverosky et al. [35] showed that M-R/S is biased towards rejecting long-range depen- 
dence in the process when the high number of lags is used. The bias has been shown on the 
behavior of the estimated rescaled ranges. In our simulations, we examine whether such a bias is 
also present for the estimation of Hurst exponent H. 

Finite sample properties of the rescaled range analysis and the modified rescaled range analysis 
were discussed in several papers. Most importantly, the condition for the process to be character- 
ized as uncorrelatcd, which is H = 0.5, was shown to hold only asymptotically. For finite samples, 
the estimated Hurst exponent can differ significantly. Anis & Lloyd [4] showed that for the finite 
samples, expected value of the rescaled range behaves as 



Further, Peters [32] proposed to adjust the expected value of the rescaled range, which was 
claimed to be better for low scales v < 50, to 



However, Couillard & Davison [11] and Grech & Mazur [15] showed that the expected values 
based on the Monte Carlo simulations are closer to the method of Anis & Lloyd [4] . Importantly, 
they also showed that the standard deviation of the estimates does not converge according to the 
central limit theorem, which is as a square root of the time series length. Such finding is essential 
for the confidence interval construction as the slower convergence of standard deviation makes the 
confidence intervals wider. 

3. Monte Carlo simulations methodology 

In this section, we briefly describe the types of simulated processes with different kinds of 
memory, depict types of distributions of the innovations £t for the simulated processes - standard 
normal, log-normal, Cauchy, log-i, gamma, inverse gamma, Laplace and log-Laplace - and specify 
the crucial parameters of Monte Carlo simulations. 

3.1. Processes 

To uncover whether R/S and M-R/S arc able to efficiently distinguish between short and long- 
range memory, we simulate three different types of processes - independent, short-range dependent 
and long-range dependent. For independent process, we simply choose i.i.d. process, i.e. Xt = £f 




(3) 




For short-range dependent process, we use an autorcgressive process with AR(1) specification, 
i.e. Xt = Oxt-i + £f And for long-range dependent process, we simulate fractionally integrated 
ARMA process ARFIMA(0, d, 0) defined as x t = a-i{d)x t -i + e t where < d < 0.5 is a free 

parameter, related to Hurst exponent as H = d + 0.5, and a,i(d) = dT(i — d)/(T(l — d)T(l + i)). 
Infinity is trimmed to 100 lags in the simulations. Note again that the noise terms et are drawn 
from various distributions described in the following subsection. Process X t is then taken as an 
integrated series -X" t = £j for i = 1, . . . , T and Xq = 0. 

3.2. Distributions description 

Processes based on normal distribution are used in a majority of models for simplicity. We use 
iV(0, 1) as a benchmark and for a comparison as well since majority of sampling properties studies 
base the results on such process. 

Log-normal distribution has several properties which are in hand with stylized facts such as the 
fat tails and the skewness of e.g. financial returns [10]. For X ~ A r (0, 1), Y = e x is log-normally 
distributed. Moreover, we subtract a unity from Y so that the distribution has zero mean and a 
minimum of —1. We apply the subtraction to all distributions to have the mean of zero. 

Cauchy distribution is a stable distribution with a = 1. For our purposes, we use the standard 
Cauchy distribution, which is equal to the Student's t distribution with one degree of freedom, 
defined on the basis of the stable distributions as 5(1, 0, y/2/2, 0). As Cauchy distribution is 
symmetric, we isolate the effect of the heavy tails from the one of the skewness. 

Log-t distribution has similar properties to log-normal distribution. However, log-i distribution 
is more skewed and has heavier tails than log-normal distribution. If X ~ t n is Student's t 
distributed with n degrees of freedom, then Y = e x is log-t distributed with n degrees of freedom. 
In simulations, we use n = 5 to have the final distribution with finite mean, variance, skewness 
and kurtosis. 

Gamma distribution is described by parameters k and 9 and can be seen as a sum of k inde- 
pendent exponential distributions with mean 9 for k € N. Such gamma distribution has a mean 
of k9, a standard deviation of \Jk9 2 , a skewness of 2/\/fe, an excess kurtosis of 6/k and is noted 
as G(k9, Vk9 2 ). For the simulations, we set k = 4 and 9 = 0.25 and subtract a unity to obtain 
G(0, 2) with the minimum of — 1, the skewness of 1 and the excess kurtosis of 3/2. 

Inverse gamma distribution is obtained from an inverse of gamma distribution. If X ~ 
G(k9,Vk9 2 ), then Y = 1/X is inverse gamma distributed with a mean 9j(k — 1), a standard 
deviation of ^J9 2 /((k — \) 2 (k — 2)), a skewness of A\fk — 2/ (k — 3) and an excess kurtosis of 
(30fc — 66)/ (fc — 3)(fc — 4). We set the parameters as for the gamma distribution and subtract unity 
to obtain skewed series with an undefined kurtosis. 

Laplace distribution is similar to normal distribution but is based on absolute deviations from a 
mean rather than squared deviations. Laplace distribution is sometimes called double-exponential 
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as it is formed from two exponential distributions which equal at and are symmetric around the 
mean. The distribution is then symmetric and leptokurtic with excess kurtosis of 3. We use the 
standard Laplace distribution with zero mean and a standard deviation of one - L(0, 1). 

Log-Laplace distribution has again similar properties to the other ones based on exponential of 
the symmetric distributions - a nonzero skewness and an excess kurtosis. If X ~ £(0, 1) is Laplace 
distributed, then Y — e x is log-Laplace distributed. As for the other exponential transformations, 
we subtract a unity to arrive at zero mean. 

3.3. Choice of parameters 

The choice of parameters is crucial for the final results of the simulations and the estimates. 
The sub- length v is set equal to the power of a fixed integer value [37]. We use a basis b = 2, 
a minimum power p m in and a maximum power p m ax so that we get v = b Pmin , . . . , b Pma * , where 
b Pmin > b is a minimum scale and b Pma:c < T is a maximum scale. We use the minimum scale 
v m in = 2 5 and the maximum scale v max = T. Such choice avoids biases of the low scales where the 
standard deviations can be estimated inefficiently [32, 15, 28]. The used intervals are contingent 
and non-overlapping (see [13] for discussion). 

In the following section, we simulate 1,000 time series with lengths ranging from 2 9 to 2 14 
for the standard normal distribution and estimate Hurst exponent based on R/S and M-R/S for 
each series. We apply the same procedure for log-normal, Cauchy, log-i, gamma, inverse gamma, 
Laplace and log-Laplace distributions. For each distribution, we construct three different types 
of processes - independent, ARFIMA with d = 0.25 and AR(1) with 9 = 0.25 - and comment 
on an accuracy of the methods for uncorrelated, short-term correlated and long-term correlated 
processes. To compare the quality of estimators for different types of processes and distribution 
of innovations, we compare their bias 8 = (Hi — E(H)), variance a 2 H = (H 2 ) — (Hi) 2 and mean 
squared error MSE = ((Hi — E(H)) 2 ) = 6 2 + a 2 Hl where Hi is the estimated Hurst exponent 
for the simulated series, E(.) stands for the expected value, and () is an average operator for 
i = 1, . . . , 1000. Note that the expected value of H for independent and short-range dependent 
processes is not 0.5 but the value needs to be calculated from Eq. 3. For long-range dependent 
processes, we expect H = 0.75 for ARFIMA(0, 0.25,0). These expected values are used to calculate 
bias and MSE in the following section. 

4. Monte Carlo simulations results 

The results of Monte Carlo simulations arc summarized in Tables 1-6. We now discuss 
the results for the three types of processes separately. For i.i.d. processes, we expect that the 
methods yield same results. For long-range dependent processes, the methods should give similar 
results. However, M-R/S is expected to underestimate the true Hurst exponent because long- 
range dependence can be misinterpreted by the method as a strong short-range dependence. Such 
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underestimation should be only small to keep the method usable. For short memory processes, 
R/S is expected to be biased upwards while M-R/S is expected to be unbiased. 

4-1. Independent process 

Tables 1 and 2 contain the results for i.i.d. simulations for R/S and M-R/S, respectively. 
Regarding the bias, there are three interesting findings. First, R/S tends to overestimate the Hurst 
exponent, while M-R/S method underestimates it. However, the biases are of an order of 10~ 3 and 
thus rather negligible. Second, the processes based on logarithmic distributions lead to stronger 
underestimation of the Hurst exponents. This is true for log-i and log-Laplace distributions but 
not for log-normal distribution for both R/S and M-R/S. As the distributions have both non-zero 
skewness and excess kurtosis, the bias is an implication of one or both of theses. Similar yet not so 
strong bias is observed for Cauchy distribution. As Cauchy distribution has infinite variance (and 
is symmetric even though it has an undefined skewness), the downward bias seems to be caused 
by a combination of excess kurtosis and strong skewness. However, the values of skewness and 
kurtosis have to be rather high to cause such a bias since other distributions (log-normal, Laplace, 
gamma and inverse-gamma) are not affected in such a way. Third, the bias in general does not 
show any tendency to decrease with an increasing time series length. 

As the biases are rather small, the qualitative results for variance and MSE of the estimators 
practically overlay. Here, we find two important outcomes. First, the estimators for log-i and 
log-Laplace distributions show the smallest variance and MSE. Second, the variance and MSE of 
the estimators decrease with an increasing time series length for all distributions and both rescaled 
range methods. Overall, the two methods do not differ significantly for i.i.d. processes analyzed 
here. 

4.2. ARFIMA fO, 0.25, 0) 

The results are more interesting for the ARFIMA processes (Tables 3 and 4). For both R/S 
and M-R/S, the bias of estimators increases (in absolute terms) with a time series length. For 
example for R/S and a standard normal distribution, the estimator is unbiased for a length of 2 9 
but is strongly biased for 2 14 with a bias of -0.06. This finding docs not change for any distribution 
used. Such a result is most likely caused by the fact that the rescaled range analysis does not 
have an expected value of 0.75 for ARFIMA process with d — 0.25. Unfortunately, these expected 
values are not yet analytically solved. As for the case of independent processes, when the expected 
values are derived from Eq. 3, the expected value of Hurst exponent decreases with increasing time 
series length. This is reflected in the increasing bias of the estimators for ARFIMA simulations. 

More importantly, the bias differs remarkably for R/S and M-R/S. For the standard normal 
distribution with the shortest length, the bias is 100 times higher for M-R/S than for R/S. Even 
though this is a rather extreme case majorly caused by the fact that R/S is almost unbiased for the 
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normal distribution, the differences arc quite large even for the other distributions. R/S estimator 
is most biased for log-t distribution, yet M-R/S estimator is more than 5 times more biased. For 
the others, the bias of M-R/S is around 20 times the bias of R/S. This huge disproportion decreases 
with an increasing time series length. Nevertheless, the M-R/S bias is between 2-3 times higher 
than the one of R/S even for the length of 2 14 for all distributions. 

Moving to the absolute terms, the bias of R/S is insignificant for the shortest time series for 
all distributions and increases to around -0.06 for the longest time series for all but the \og-t 
distribution with a bias of -0.0969. For M-R/S, the bias is smallest for the shortest series (around 
-0.1) and increases to approximately -0.15 for the series of 2 observations for all examined 
distributions. This is a very important finding which indicates that M-R/S misinterprets the long- 
memory of the process as the short-memory. Even though the method is constructed to rid of the 
strong short-term correlations, our results show that the correction for the short memory is way 
too strong. This practically disqualifies the method from being used for long memory analysis. 

The variances of the M-R/S estimator are lower than the ones of the R/S for all distributions 
and all time series lengths (the R/S variances are around 0.5 times higher than these of M-R/S). 
The variance decreases with an increasing series length for both methods and all distributions. 
Again, the variances for log-i and log-Laplace distributions are somewhat lower. Even though the 
variances of M-R/S are lower, the effect of the bias is so strong that the MSE of M-R/S is higher 
than the one of R/S (approximately 0.8 times higher). The difference is the lowest for the log-t 
distribution. Interestingly, the bias increases the MSE so much that the MSE increases with a time 
series length for M-R/S. For R/S, the MSE standardly decreases with the number of observations. 

4.3. AR(1) 

For short-range dependent AR(1), the results are quite as expected but not as strong. Even 
though R/S overestimates H for a strong majority of distributions and lengths of the series, the 
maximum bias is 0.0309 for Laplace distribution of length 2 9 . The bias decreases with the number 
of observations for all but log-t distribution. This is caused by the fact that the estimator is 
unbiased for log-t distribution for small samples. For M-R/S, the estimates are less biased than 
the ones of R/S for majority of distributions. However, M-R/S is more biased for log-t and log- 
Laplacc distributions for all sample lengths. Moreover, the bias of M-R/S does not decrease with 
sample length in all cases. Similarly to the i.i.d. case, the estimates are pushed downwards for 
\og-t and log-Laplace distributions. While this is true for both R/S and M-R/S, the effect is 
different. For R/S, the bias is pushed from the positive numbers closer to zero. For M-R/S, the 
bias is negative and pushed further away from zero. 

Even though MSE is dominated by variance for M-R/S, the effect of the bias is stronger for 
R/S. Nevertheless, the overall qualitative results remain the same for both the variance and MSE 
of the estimators. The variance and MSE of M-R/S are lower for shorter series. However, the 
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difference becomes negligible for longer series. Note that again the variance and MSE are lower 
for log-i and log-Laplace distributions. 

5. Conclusions 

In this paper, we analyzed the robustness of rescaled range analysis and modified rescaled 
range analysis to different distributional and memory properties. The bias, variance and mean 
squared errors of the estimators have been studied for independent, short-range dependent and 
long-range dependent processes with innovations from 8 different distributions. The simulations 
study uncovered mainly four interesting implications. 

First, extreme excess kurtosis combined with high skewness cause the estimates to be downward 
biased. Second, M-R/S in its simplest way is not able to distinguish between short and long 
memory efficiently. This is a very crucial result. It implies that M-R/S should not be used for 
Hurst exponent estimation in presence of short-term memory in the underlying process. This 
result is in hand with Tevcrosky et al. [35], who showed that M-R/S is not suitable for adjusting 
the rescaled ranges. However, this should not disqualify the method altogether. It shows that the 
method needs to be used properly and in conjunction with other methods. One of the examples 
is using M-R/S with bootstrapped confidence intervals as has been already done several times 
[19, 31]. This way, both types of dependence can be efficiently distinguished. Third, R/S is biased 
by the short-range dependence in the series but in a lesser way than standardly claimed. On 
contrary, R/S is less biased than M-R/S for some distributions. And fourth, apart from very 
extreme cases of high leptokurtosis and skewness, both methods are quite robust to different 
distributional properties. 
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Table 1: Results of Monte Carlo simulations for R/S — H=0.5 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


0.0034 


0.0007 


0.0038 


0.0011 


0.0005 


0.0015 




log-normal 


0.0003 


-0.0026 


0.0014 


0.0038 


0.0042 


0.0053 




log-t 


-0.0262 


-0.0219 


-0.0208 


-0.0209 


-0.0206 


-0.0200 




Cauchy 


-0.0088 


-0.0090 


-0.0049 


0.0012 


0.0028 


0.0052 




Laplace 


0.0043 


0.0019 


0.0069 


0.0050 


0.0060 


0.0050 




log-Laplace 


-0.0154 


-0.0118 


-0.0110 


-0.0099 


-0.0065 


-0.0066 




gamma 


0.0000 


0.0028 


0.0055 


0.0017 


0.0009 


0.0007 




inv-gamma 


-0.0002 


0.0023 


0.0019 


0.0041 


0.0047 


0.0030 


variance 


normal 


0.0070 


0.0036 


0.0024 


0.0015 


0.0010 


0.0008 




log-normal 


0.0062 


0.0036 


0.0022 


0.0016 


0.0011 


0.0007 




log-t 


0.0033 


0.0019 


0.0011 


0.0006 


0.0004 


0.0002 




Cauchy 


0.0044 


0.0027 


0.0018 


0.0013 


0.0009 


0.0006 




Laplace 


0.0063 


0.0038 


0.0023 


0.0016 


0.0010 


0.0007 




log-Laplace 


0.0044 


0.0027 


0.0017 


0.0011 


0.0008 


0.0006 




gamma 


0.0058 


0.0036 


0.0023 


0.0014 


0.0011 


0.0008 




inv-gamma 


0.0058 


0.0035 


0.0022 


0.0014 


0.0011 


0.0008 


MSE 


normal 


0.0070 


0.0036 


0.0024 


0.0015 


0.0010 


0.0008 




log-normal 


0.0061 


0.0036 


0.0022 


0.0016 


0.0011 


0.0007 




log-t 


0.0040 


0.0024 


0.0015 


0.0011 


0.0008 


0.0006 




Cauchy 


0.0045 


0.0028 


0.0018 


0.0013 


0.0010 


0.0007 




Laplace 


0.0064 


0.0038 


0.0023 


0.0017 


0.0011 


0.0008 




log-Laplace 


0.0046 


0.0028 


0.0018 


0.0012 


0.0008 


0.0006 




gamma 


0.0058 


0.0036 


0.0024 


0.0014 


0.0011 


0.0008 




inv-gamma 


0.0058 


0.0035 


0.0022 


0.0014 


0.0011 


0.0008 



Table 2: Results of Monte Carlo simulations for M-R/S - H=0.5 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


-0.0071 


-0.0068 


-0.0017 


-0.0032 


-0.0030 


-0.0015 




log-normal 


-0.0085 


-0.0084 


-0.0030 


0.0005 


0.0016 


0.0031 




log-t 


-0.0336 


-0.0266 


-0.0241 


-0.0234 


-0.0225 


-0.0215 




Cauchy 


-0.0183 


-0.0160 


-0.0100 


-0.0029 


-0.0005 


0.0022 




Laplace 


-0.0070 


-0.0059 


0.0010 


0.0006 


0.0024 


0.0020 




log-Laplace 


-0.0235 


-0.0171 


-0.0149 


-0.0129 


-0.0088 


-0.0084 




gamma 


-0.0105 


-0.0045 


0.0001 


-0.0025 


-0.0025 


-0.0022 




inv-gamma 


-0.0098 


-0.0046 


-0.0028 


0.0003 


0.0017 


0.0005 


variance 


normal 


0.0067 


0.0035 


0.0024 


0.0015 


0.0010 


0.0008 




log-normal 


0.0060 


0.0035 


0.0021 


0.0015 


0.0011 


0.0007 




log-t 


0.0034 


0.0018 


0.0011 


0.0006 


0.0004 


0.0002 




Cauchy 


0.0043 


0.0027 


0.0018 


0.0013 


0.0009 


0.0006 




Laplace 


0.0064 


0.0038 


0.0023 


0.0016 


0.0010 


0.0007 




log-Laplace 


0.0043 


0.0027 


0.0017 


0.0011 


0.0008 


0.0006 




gamma 


0.0057 


0.0035 


0.0023 


0.0014 


0.0011 


0.0008 




inv-gamma 


0.0056 


0.0034 


0.0022 


0.0014 


0.0011 


0.0007 


MSE 


normal 


0.0068 


0.0036 


0.0024 


0.0015 


0.0010 


0.0008 




log-normal 


0.0061 


0.0036 


0.0021 


0.0015 


0.0011 


0.0007 




log-t 


0.0046 


0.0025 


0.0016 


0.0012 


0.0009 


0.0007 




Cauchy 


0.0047 


0.0030 


0.0019 


0.0013 


0.0009 


0.0006 




Laplace 


0.0064 


0.0038 


0.0023 


0.0016 


0.0010 


0.0007 




log-Laplace 


0.0049 


0.0030 


0.0019 


0.0013 


0.0008 


0.0006 




gamma 


0.0058 


0.0035 


0.0023 


0.0014 


0.0011 


0.0008 




inv-gamma 


0.0057 


0.0034 


0.0022 


0.0014 


0.0011 


0.0007 
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Table 3: Results of Monte Carlo simulations for R/S - H=0.75 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


-0.0010 


-0.0031 


-0.0124 


-0.0278 


-0.0439 


-0.0600 




log-normal 


-0.0021 


-0.0005 


-0.0121 


-0.0266 


-0.0433 


-0.0618 




log-t 


-0.0170 


-0.0287 


-0.0384 


-0.0603 


-0.0786 


-0.0969 




Cauchy 


-0.0041 


-0.0022 


-0.0189 


-0.0344 


-0.0512 


-0.0656 




Laplace 


0.0016 


0.0013 


-0.0127 


-0.0262 


-0.0437 


-0.0608 




log-Laplace 


-0.0041 


-0.0119 


-0.0208 


-0.0421 


-0.0569 


-0.0745 




gamma 


0.0029 


-0.0001 


-0.0162 


-0.0282 


-0.0449 


-0.0596 




inv-gamma 


-0.0050 


-0.0045 


-0.0102 


-0.0272 


-0.0445 


-0.0606 


variance 


normal 


0.0109 


0.0057 


0.0032 


0.0021 


0.0013 


0.0008 




log-normal 


0.0096 


0.0054 


0.0034 


0.0020 


0.0013 


0.0008 




log-t 


0.0089 


0.0041 


0.0022 


0.0011 


0.0007 


0.0004 




Cauchy 


0.0079 


0.0039 


0.0025 


0.0015 


0.0011 


0.0008 




Laplace 


0.0102 


0.0057 


0.0033 


0.0022 


0.0013 


0.0009 




log-Laplace 


0.0083 


0.0046 


0.0027 


0.0015 


0.0009 


0.0006 




gamma 


0.0093 


0.0054 


0.0032 


0.0022 


0.0013 


0.0009 




inv-gamma 


0.0099 


0.0052 


0.0034 


0.0021 


0.0014 


0.0008 


MSE 


normal 


0.0109 


0.0057 


0.0034 


0.0029 


0.0032 


0.0044 




log-normal 


0.0096 


0.0054 


0.0035 


0.0027 


0.0032 


0.0046 




log-t 


0.0092 


0.0049 


0.0037 


0.0048 


0.0068 


0.0097 




Cauchy 


0.0079 


0.0039 


0.0028 


0.0027 


0.0037 


0.0051 




Laplace 


0.0102 


0.0057 


0.0035 


0.0029 


0.0032 


0.0046 




log-Laplace 


0.0083 


0.0047 


0.0031 


0.0033 


0.0041 


0.0062 




gamma 


0.0093 


0.0054 


0.0035 


0.0029 


0.0033 


0.0044 




inv-gamma 


0.0099 


0.0052 


0.0035 


0.0028 


0.0033 


0.0045 



Table 4: Results of Monte Carlo simulations for M-R/S - H=0.75 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


-0.1004 


-0.1026 


-0.1107 


-0.1249 


-0.1390 


-0.1504 




log-normal 


-0.0994 


-0.1006 


-0.1098 


-0.1224 


-0.1374 


-0.1513 




log-t 


-0.0896 


-0.0947 


-0.1072 


-0.1258 


-0.1446 


-0.1605 




Cauchy 


-0.0970 


-0.0991 


-0.1124 


-0.1272 


-0.1421 


-0.1528 




Laplace 


-0.0982 


-0.0988 


-0.1106 


-0.1230 


-0.1386 


-0.1514 




log-Laplace 


-0.0942 


-0.1018 


-0.1103 


-0.1297 


-0.1433 


-0.1576 




gamma 


-0.0963 


-0.1000 


-0.1144 


-0.1251 


-0.1394 


-0.1500 




inv-gamma 


-0.1020 


-0.1034 


-0.1089 


-0.1239 


-0.1386 


-0.1506 


variance 


normal 


0.0074 


0.0040 


0.0024 


0.0017 


0.0011 


0.0007 




log-normal 


0.0067 


0.0038 


0.0025 


0.0016 


0.0011 


0.0007 




log-t 


0.0042 


0.0021 


0.0012 


0.0006 


0.0003 


0.0002 




Cauchy 


0.0051 


0.0027 


0.0019 


0.0013 


0.0009 


0.0007 




Laplace 


0.0069 


0.0040 


0.0025 


0.0018 


0.0011 


0.0008 




log-Laplace 


0.0052 


0.0030 


0.0018 


0.0012 


0.0007 


0.0005 




gamma 


0.0063 


0.0037 


0.0025 


0.0017 


0.0011 


0.0008 




inv-gamma 


0.0068 


0.0037 


0.0025 


0.0017 


0.0011 


0.0007 


MSE 


normal 


0.0175 


0.0145 


0.0146 


0.0173 


0.0204 


0.0234 




log-normal 


0.0165 


0.0139 


0.0146 


0.0166 


0.0200 


0.0236 




log-t 


0.0123 


0.0110 


0.0127 


0.0164 


0.0212 


0.0260 




Cauchy 


0.0146 


0.0125 


0.0145 


0.0174 


0.0211 


0.0240 




Laplace 


0.0165 


0.0138 


0.0147 


0.0169 


0.0203 


0.0237 




log-Laplace 


0.0141 


0.0134 


0.0140 


0.0180 


0.0212 


0.0254 




gamma 


0.0156 


0.0137 


0.0155 


0.0174 


0.0206 


0.0233 




inv-gamma 


0.0172 


0.0144 


0.0144 


0.0170 


0.0204 


0.0234 
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Table 5: Results of Monte Carlo simulations for R/S - AR(1) 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


0.0288 


0.0220 


0.0216 


0.0192 


0.0173 


0.0134 




log-normal 


0.0253 


0.0208 


0.0180 


0.0205 


0.0194 


0.0158 




log-* 


0.0042 


0.0022 


-0.0026 


-0.0048 


-0.0055 


-0.0063 




Cauchy 


0.0158 


0.0120 


0.0139 


0.0145 


0.0165 


0.0157 




Laplace 


0.0309 


0.0256 


0.0232 


0.0209 


0.0186 


0.0165 




log-Laplace 


0.0123 


0.0094 


0.0072 


0.0090 


0.0058 


0.0059 




gamma 


0.0253 


0.0220 


0.0204 


0.0177 


0.0146 


0.0149 




inv-gamma 


0.0273 


0.0226 


0.0222 


0.0186 


0.0175 


0.0164 


variance 


normal 


0.0066 


0.0040 


0.0023 


0.0015 


0.0010 


0.0007 




log-normal 


0.0061 


0.0039 


0.0024 


0.0015 


0.0011 


0.0008 




Iog-t 


0.0036 


0.0018 


0.0010 


0.0005 


0.0004 


0.0002 




Cauchy 


0.0045 


0.0029 


0.0018 


0.0014 


0.0009 


0.0007 




Laplace 


0.0064 


0.0040 


0.0025 


0.0015 


0.0011 


0.0009 




log-Laplace 


0.0044 


0.0026 


0.0017 


0.0010 


0.0007 


0.0006 




gamma 


0.0066 


0.0037 


0.0023 


0.0016 


0.0011 


0.0008 




inv-gamma 


0.0060 


0.0038 


0.0024 


0.0017 


0.0010 


0.0008 


MSE 


normal 


0.0074 


0.0045 


0.0028 


0.0019 


0.0013 


0.0009 




log-normal 


0.0067 


0.0043 


0.0027 


0.0019 


0.0014 


0.0010 




Iog-t 


0.0036 


0.0018 


0.0010 


0.0005 


0.0004 


0.0003 




Cauchy 


0.0047 


0.0030 


0.0020 


0.0016 


0.0011 


0.0009 




Laplace 


0.0073 


0.0046 


0.0030 


0.0019 


0.0014 


0.0011 




log-Laplace 


0.0045 


0.0027 


0.0018 


0.0011 


0.0007 


0.0006 




gamma 


0.0072 


0.0042 


0.0027 


0.0019 


0.0013 


0.0010 




inv-gamma 


0.0068 


0.0043 


0.0028 


0.0020 


0.0013 


0.0010 



Table 6: Results of Monte Carlo simulations for M-R/S - AR(1) 





512 


1024 


2048 


4096 


8192 


16384 


bias 


normal 


-0.0019 


-0.0052 


-0.0019 


-0.0023 


-0.0022 


-0.0044 




log-normal 


-0.0055 


-0.0066 


-0.0061 


-0.0015 


-0.0004 


-0.0019 




log-* 


-0.0246 


-0.0244 


-0.0256 


-0.0255 


-0.0246 


-0.0235 




Cauchy 


-0.0143 


-0.0139 


-0.0095 


-0.0071 


-0.0030 


-0.0022 




Laplace 


0.0006 


-0.0011 


-0.0009 


-0.0011 


-0.0009 


-0.0014 




log-Laplace 


-0.0182 


-0.0170 


-0.0170 


-0.0124 


-0.0135 


-0.0115 




gamma 


-0.0056 


-0.0050 


-0.0037 


-0.0044 


-0.0049 


-0.0029 




inv-gamma 


-0.0034 


-0.0049 


-0.0016 


-0.0032 


-0.0020 


-0.0016 


variance 


normal 


0.0060 


0.0037 


0.0022 


0.0015 


0.0010 


0.0007 




log-normal 


0.0055 


0.0036 


0.0023 


0.0014 


0.0010 


0.0008 




log-* 


0.0033 


0.0016 


0.0009 


0.0005 


0.0004 


0.0002 




Cauchy 


0.0041 


0.0027 


0.0017 


0.0013 


0.0008 


0.0007 




Laplace 


0.0059 


0.0037 


0.0023 


0.0014 


0.0010 


0.0008 




log-Laplace 


0.0040 


0.0024 


0.0016 


0.0010 


0.0007 


0.0006 




gamma 


0.0060 


0.0035 


0.0022 


0.0016 


0.0011 


0.0007 




inv-gamma 


0.0054 


0.0035 


0.0022 


0.0016 


0.0010 


0.0008 


MSE 


normal 


0.0060 


0.0037 


0.0022 


0.0015 


0.0010 


0.0007 




log-normal 


0.0055 


0.0037 


0.0023 


0.0014 


0.0010 


0.0008 




log-t 


0.0039 


0.0022 


0.0016 


0.0012 


0.0010 


0.0008 




Cauchy 


0.0043 


0.0029 


0.0018 


0.0014 


0.0008 


0.0007 




Laplace 


0.0059 


0.0037 


0.0023 


0.0014 


0.0010 


0.0008 




log-Laplace 


0.0043 


0.0027 


0.0019 


0.0011 


0.0009 


0.0007 




gamma 


0.0060 


0.0035 


0.0022 


0.0016 


0.0011 


0.0007 




inv-gamma 


0.0054 


0.0035 


0.0022 


0.0016 


0.0010 


0.0008 
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